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We derive spectral representations for the different components of the 4-point function at finite 
temperature in the real time formalism in terms of five real spectral densities. We explicitly calculate 
all these functions in QED in the hard thermal loop approximation. The Ward identities obeyed by 
the 1-loop 3- and 4-point functions in real time and their spectral functions are derived. We compare 
■ our results with those derived previously in the imaginary-time formalism for retarded functions in 

hot QCD, and we discuss the generalization of our results to non-equilibrium situations. 
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I. INTRODUCTION 



One of the most significant advances in quantum field theory at finite temperature is the effective perturbation 
theory of Braaten and Pisarski which is based on the resummation of hard thermal loops (HTLs) . HTLs were 
originally obtained by computing one-loop diagrams in the imaginary time formulation of finite temperature field 
theory. HTLs are ultraviolet finite, gauge invariant and satisfy simple Ward identities fp rfl- These remarkable 
properties have been traced to the fact that HTLs describe simple semiclassical physics In particular, the 

HTLs can be derived from a set of classical kinetic equations [|||l0| as well as from a non-local effective Lagrangian 
' |§- The computation of HTLs is fairly technical because of their non-trivial momentum and energy dependence, but 
can be simplified by using the Ward identities 

The analytical structure of Green functions at zero or non-zero temperature, in the imaginary or real-time formula- 
fS |' tions, is controlled by their spectral densities, combined with asymptotic boundary conditions. At finite temperature 
the spectral functions can also be directly related to transport coefficients p^ , |l3[ . The study of spectral functions 
helps us to understand the quasi-particle structure of a field theory and to identify the microscopic processes underly- 
ing the dynamics. But Green functions and their spectral functions are not easy quantities to evaluate perturbatively 
L* ' at non-zero temperature, especially for many-point functions. The two-point functions and their spectral densities 
have been widely studied and applied to quark gluon plasma (QGP) investigations [1^,^,^,^2|. Our knowledge about 
■ many-point spectral densities is still much less comprehensive p~5|— P0|| . T he spectral 



representations of 3-point func- 
tions for self- interacting scalar fields were discussed in Refs. [|15 ,|16|,|19 20 1. In |18j 3- and 4-point spectral densities for 
pure gluon dynamics were calculated in the HTL approximation using the imaginary time formalism (ITF) p|, and 
in the same approximation the 3-point spectral density in QED has been calculated within the real time formalism 
(RTF) §}. " _ " 

The popularity of the RTF is diminished by the technical complications resulting from the doubling of degrees of 
freedom [ pl"| . However, in the ITF, one has to perform an analytic continuation of the imaginary external energy 
variables to the real axis, which is avoided in the RTF. In addition, the ITF is restricted to equilibrium situations, 
while the RTF can be extended to investigate non-equilibrium systems J2l"|-p3| . 

In this paper we study 4-point functions in QED using the real-time formalism. We adopt the Keldysh or closed 
time path (CTP) contour p2|. In SecJO we introduce the Keldysh representation. The spectral representations of 



the 4-point functions is derived in Sec. ITJ. In Sec. IV we use the Keldysh representation to calculate the 4-point HTL 
vertex functions in hot QED. The 4-point spectral densities in the HTL approximation are extracted in Sec. [V and 
shown to degenerate to a single real spectral function. In Sec. [vi| we derive real-time Ward identities between the 3- 
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and 4-point HTLs and their s pect ral functio ns. T he generalization of our results to certain types of non-equilibrium 
situations is discussed in Sec. VII, while Sec. VIII gives a summary of our results. 



II. KELDYSH REPRESENTATION 



We will use the Keldysh representation [£l|,^2| of the real-time formalism (RTF) throughout this paper. In this 
representation the single-particle propagator for free bosons has in momentum space the form [^l|p2[ 

D{K) = ( T^^i °, ) 2mS{K 2 -rn 2 ) ( ... ^ k °\ h , + n f ^ ) (1) 

V K*-rn»-ie J J \ 6{k ) + n B {k ) n B (k Q ) J 

where K — (ka,k), k — |fe|, 9 denotes the step function, and the equilibrium distribution function is given by 
n s(fco) = l/[exp(|fco|/T) — !]• For fermions the bare propagator can be written as 



S(K) = {If + m) 



V K»-m*-U J J \0{k Q ) -n F {k ) n F (k ) 



(2) 



where the Fermi-Dirac distribution is given by np(ko) = l/[exp(|fco|/T) + 1]. The components of these propagators 
are not independent, but fulfill the relation 

Gn — G12 — G21 + G22 = 0, (3) 

where G stands for D or S, respectively. 

By an orthogonal transformation of these 2x2 matrices one arrives at the representation of the propagators in 
terms of advanced and retarded propagators which was first introduced by Keldysh p3). The three nonvanishing 
components of this representation are |21[] 

Gr = Gn — G12, 
Ga = Gn — G21; 

G F = Gn + G*22- (4) 



The inverted relations read 



Gu = 2 (Gf + Ga + Gr), 
G12 — — {Gf + Ga — Gr), 
G21 — 2 (Gf — Ga + Gr), 

G22 — 2 {Gf - Ga - Gr). (5) 



The relations (0-0) also hold for full propagators. 

Using ([j]) and (g) in ([|), the free propagators are given in the Keldysh representation by 

Dr{K) = ' 



K 2 — m 2 + i sgn{ko)e 

D A {K) = — - ; — — , 

K A — m l — zsgn(Ko)e 

D F {K) = (l+2nfl(*o)) sgn(fc ) (p R {K) - D A (K)) = -2m (l+2n B (fco)) 5{K 2 -m 2 ) (6) 



for bosons and 



Sr{K) = 




If + m 




K 2 


— m 2 + i sgn(fco)e 


S A (K) = 




If + m 


K 2 


— m 2 — i sgn(fco)e 


S F {K) = 


(1- 


2n F {k )) sgn(fco) ( 



{If + m)D R {K) , 
{If + m)D A {K), 



= -2ni {If + m) (l-2n F (k )) S{K 2 - m 2 ) = {If + m)D F {K) (7) 
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for fermions, where the last equation in each of (|6|) and (|7|) is a consequence of the dissipation-fluctuation theorem. 
Although Da.r — Da,r we keep the bar in our notation to more easily recognize fermion and boson propagators in 
the calculations below. 



III. SPECTRAL REPRESENTATION OF THE 4-POINT VERTEX 



In this section we briefly review some useful relations among the different thermal components of the 4-point 
functions, and we derive their spectral representations. Some of these relations have already been reported in the 
literature @-||,||| using different notation. We discuss these results at the end of this section. For any 1PI four- 
point function we define the four incoming external momenta as Pi, P2, P3, and P4 = — (P1+P2+P3). In this section 
we suppress all indices except Keldysh indices. 



A. Largest and smallest time equations 

If ti is the largest time argument, one can obtain the "largest time equations" [ fi^]l5p4| ] : 

#43 #32 #21 (Gabcl + Gabc2) = #43 #31 #12 (Gabcl + G a bc2) = , (8a 

#42 #23 #31 (Gabcl + G a bc2) = #42 #23 #13 (G a bcl + G a bc2) = , (8b 

#41 #13 #32 (Gabcl + G a bc2) = #41 #12 #23 (G a bcl + G abc2) = , (8c 

#32 #21 #14 (Gable + G a b2c) = #31 #12 #24 (G able + G a b2c) = , (8d 

#23 #34 #41 (Galbc + G a 2bc) = #24 #43 #31 (G a lbc + G a 2bc) = , (8e 

#13 #32 #24 (Gi Q b c + G2abc) = #12 #23 #34 (G\ a bc + G2abc) = , (8f 

where a, b, c can be either 1 or 2, and 9ij = 9(ti — tj) is the step function. 

By tilde conjugation (which replaces time ordering with anti-time ordering pl| , p3[ ) one obtains equations of the 
form: 

#12 #23 #34 (G aoc i + G a b c 2) = #13 #32 #24 (G a bcl + G a bc2) = , (9a) 

#41 #12 #23 (Gable + G a b2c) = #42 #21 #13 (G a blc + G a b2c) = , (9b) 

#14 #43 #32 (Galbc + G a 2bc) = #13 #34 #42 (G a ib c + G a 2bc) = , (9c) 

#42 #23 #31 (Glabc + G2abc) = #43 #32 #21 (Gi a bc + G2abc) = , (9d) 

Eqs. (|h and @ are the analogues of the "largest time equations" and "smallest time equations", respectively, of 

Rcf . J14| . They will be used extensively in the derivation of the spectral representations below. Their generalization 
to arbitrary n-point functions is straightforward. 



B. Derivation of spectral representations 



One can construct the "retarded-advanced" vertex functions from the sixteen components of the real-time 4-point 
function as in [^9 21 1. The KMS conditions allows to reduce them to the following seven combinations [ p5p^ ]: 

(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 
(lOf) 
(10g) 

It is straightforward to show that the first four are the standard retarded vertex functions 

G r aaa = G m ] G araa = Gfl2 j G a ara = Gfl3 ; G aaar = G R4 . (11) 



Graaa 


= Gini - 


- G1211 - 


- G1121 - 


- G1112 - 


- G1122 - 


- G1221 - 


- G1212 - 


- G1222 


Garaa 


= Gini - 


- G2111 - 


- G1121 - 


- G1112 - 


- G1122 - 


- G2121 - 


- G2112 - 


- G2122 


Gaara 


= Gini - 


- G2111 - 


- G1211 - 


- G1112 - 


- G1212 - 


- G2211 - 


- G2112 - 


- G2212 


Gaaar 


= Gini - 


- G2111 - 


- G1211 - 


- G1121 - 


- G1221 - 


- G2211 - 


- G2121 - 


- G2221 


Grara 


= Gini - 


- G1211 - 


- G1112 - 


- G1212 - 


- G2121 - 


- G2122 - 


- G2221 - 


- G2222 


Grraa 


= Gini - 


- G1112 - 


- G1121 - 


- G1122 - 


- G2211 - 


- G2212 - 


- G2221 - 


- G2222 


Graar 


= Gini - 


- G1121 - 


- G1211 - 


- G1221 - 


- G2112 - 


- G2122 - 


- G2212 - 


- G2222 
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The last three we will call "mixed retarded-advanced" functions. 

There are eight other vertices which one can obtain from ([Ic]) using the KMS conditions f25| , ^6f . We write down 
only the three which we need: 

(12a) 
(12b) 
(12c) 

We can re-express these vertices in terms of spectral functions by making use of various properties of theta functions. 
We will use 



Garar 


= Gini - 


h G1121 - 


h G2111 - 


h G2121 " 


h G1212 " 


\- G1222 " 


h G2212 " 


h G2222 


Garra 


= Gini - 


h G1112 - 


h G2111 - 


h G2112 " 


\- G1221 " 


\- G1222 " 


1- G2221 " 


h G2222 


Gaarr 


= Gini - 


h G1211 - 


h G2111 - 


1- G2211 - 


I- G1122 - 


1- G1222 ~ 


h G2122 - 


h G2222 



23C34C24 



012 + 021 — 1, 012021 = 0, 

(i) We begin with Gri. We rewrite G^i as 

Gri — (023 + ^32) (#34 + #43) (#24 + #42)G_R1 



23 1734 • 



Inserting the identities 



into ( |l4| ) we obtain 



Using the identities 



^i3Gfli = OuGri — OuGri = Gri 



Gri — #i2#23#34G-Ri 
+ #14#43#32Gi?i 



?12#24#43G_R1 
?13#34#42Gfli 



O14842O23GR1 
813832824GR1. 



9(lijk)G R2 
0(U]k)G R3 
6(lijk)G R4 



8n8ij8jkGR 2 
8u8ij8jkGR3 
8n8ij8jkGR4 



(i,j,k = 2,3,4) 



which result from conflicting ^-functions, we have 



G RX = [0(1234) + 0(1243) 
= [0(1234) + 0(1243) 



(1324) + 0(1342) + 0(1423) + 0(1432)] {Gri - G R2 + G R3 - G R4 ) 
(1324) + 0(1342) + 0(1423) + 0(1432)] p 12 



with 



Pl2 — Gri — Gr2 + Gr3 — Gr4 

= G1211 — G2122 + 2(Gi212 — G2121) + G1112 



G 2 



Gi 



G211 



G2 



(ii) For Gr2,Grz and Gr4 one proceeds similarly and obtains 



Gr 2 = [0(2134) + 0(2143) + 0(2314) + 0(2341) + 0(2413) + 0(2431)] (-pi 2 ), 
G R3 = [0(3124) + 0(3142) + 0(3214) + 0(3241) + 0(3412) + 0(3421)] p 12 , 
G R 4 = [0(4123) + 0(4132) + 0(4213) + 0(4231) + 0(4312) + 0(4321)] (-p u ). 

(iii) The procedure for the mixed retarded-advanced four-point functions is similar. We write 

Grraa = (012 + 02l)(014 + 04l)(023 + 032)(034 + 043)(#24 + 042)G rraa 

and use 

8(3ijk)G rraa = 0, i,j,k = 1,2,4, 
8(4ijk)G rraa = 0, i,j,fe = l,3,4, 



(13) 
(14) 
(15) 

(16) 



(17) 

(18) 
(19) 



(20a) 
(20b) 
(20c) 



(21) 



(22) 



to obtain 
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Grraa = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(2134) + 0(2143) + 0(2314) + 0(2341) + 0(2413) + 0(2431)] G r 



Using the identities 



allows us to write 



6(lijk)G rraa = 0, i,j,k = 2,3,4, 
6(2ijk)G rraa = 0, i,j,k = 1,3,4 



Grraa = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(2134) + 0(2143) + 0(2314) + 0(2341) + 0(2413) + 0(2431)] p 3 

with 

P3 Grraa G T raa • 

Analogously, we can derive 

Grara = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(3124) + 0(3142) + 0(3214) + 0(3241) + 0(3412) + 0(3421)] Pi 

Graar = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(4123) + 0(4132) + 0(4213) + 0(4231) + 0(4321) + 0(4312)] p 5 



where 



PA — G ra ra G 
PS Graar G r ac r 



rani. 



Using the Fourier integral representation of the function, 

g— iQ(ti — tj) 



i r 

~2W_ ( 



Q + ie ' 

it is straightforward to derive the spectral integral representations in momentum space: 



Gri(wi,W2,W3,o;4) = j> a~l [a^ 2 (a 3 + a7) + a^ 3 (a 2 + a7) + a^(a 2 + aig )] p X2 , 
Gk 2 (wi,w 2 ,W3,w 4 ) = <j> a^a^^ag +a7) + a 23 (ar +a7) + aj 4 (a7 + a7)](-pi 2 ) 

Gfl 3 (wi,^2,^3,^4) = <j> ^3 [ a 3l( a 2 + a 4 ) + a 32( a r + a 4 ) a 34( a r + «2~)] Pl2 , 
Gij4(wi,W 2 ,W3,W 4 ) = ^ CI4 [a4i(d2 + a 3~) + a 42( a r + a 3~) + + a 2^)](^Pl2) 

G rroo (wi,w 2 ,w 3 ,a;4) = ^ (a^ [a^ (a7 + aj) + a^ 3 {a 2 + aj) + a 1 l 4 (a7 + a7)] 

+a+[a+(a7 + a7 ) + a^ 3 (a^ + aj) + a+(a7 + a7 )]) p 3 , 
G rara (vi,u}2,U3,Ui) = <j> (a^[a^ 2 (a 3 + a7) + a^a, + aj) + af 4 (ao + a 3 )] 

+«3 [«31 ( a 2 + a 4 ) + a 32 0r + a 4 ) + + a 2~)]) Pi 

Graar(uJl,UJ2,UJ3,^4:) = <j> (af [CL^ 2 (tig + d7 ) + «i3( a 7 + a 4 ) + a 14( a 2 + a 3~)] 

+aj[ai(a7 + ) + a| 2 (a7 + ajf ) + a| 3 (a7 + a 2 )]) P5 



where 



(27T) 



dnidQ 2 dft 3 (31a) 



a ? = - ( 31b ) 
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a S = TT o 7T^-> *,i = 1,2,3,4. (31c) 

The frequency arguments of the spectral functions under the integrals are pi(£li, ^2, ^3, ^4) with ^1+^2+^3+04=0. 
The spatial momenta Pi,P2,P3, and j>4=— (P1+P2+P3) are the same on both sides of these equations and have 
therefore been suppressed. 

We have written the seven vertex functions Jl(i| ) in terms of one complex spectral density (p\ 2 = pi + ip 2 ) and 
three purely imaginary spectral densities {pz, Pi, Pb)- This means that all sixteen 4-point functions can be expressed 
in terms of five real spectral densities: 

pi = Re [n^ 1 (p2o)G2i22 + 2n^ (P20 + P4o)G 2 i2i + n^ 1 (p 4 o)G 2221 - n B 1 (pio)Gi222 - n B 1 (p 30 )G 2 2i2] , (32a) 
p 2 = -Im[rt F 1 (p 2 o)G 2 i22 + ^n F 1 (p 20 + P4o)G 2 i2i + ^f 1 (P4o)G 2 221 - rip 1 (p w )Gi222 - rip 1 (p 3Q )G 2 2i2} , (32b) 

PZ = -iP3 = 2 Im G rr aa , (32c) 
P4 = -«P4 = 2 Im G raro , (32d) 

P5 = ~iP5 = 2 Im G raar , (32e) 
where we used the following relations in momentum space resulting from the KMS condition |2l| ]: 

Giaii = e te ° G* 2l22 , G 1212 = e ^™+^ G* 2121 , G m2 = e?»» G£ 221 , (33a) 
G2111 = e^ 10 G 1222 , G 1121 = eft*» G^ 212 , (33b) 

Grraa — G rraa , G r ara — G rara , G r aar — G raar . (33c) 

These results should be compared with the expressions derived in |l5f| in which the spectral representations for the 
four retarded functions Gm, i = 1, 2, 3, 4 were given in terms of three real spectral densities. A relation between these 
three spectral densities exists which can be employed to reduce the number of independent spectral functions to two; 
this is consistent with our result. Spectral representations for the mixed retarded-advanced 4-point functions have 
also been previously discussed: in JT^ ], using a different definition of the mixed retarded-advanced 4-point functions, 
the spectral representation for one of these vertices was given in terms of six spectral densities. 

We finally note that in deriving the spectral representations (Q) no use was made of the KMS condition. In this 
form, (with the spectral densities defined by (|l6|), ( p6|) and (|28|)), Eqs. (|30| ) remain valid out of thermal equilibrium. 



IV. 4-POINT VERTEX FUNCTIONS IN QED IN THE HTL APPROXIMATION 



In this section, we calculate the seven 4-point vertex functions ( |10| ) for QED in the HTL approximation. In the 
HTL approximation QED has two non-zero 4-point vertex functions, shown in Fig. [1] below: 




G 




The other two 4-point vertex functions (four-photon vertex and four-fermion vertex) are zero in the HTL approximation 
H . We calculate the diagrams in Fig. [1] as follows: the real-time Green functions are written in the Keldysh formalism 
using (|) and @. The photon propagator in the Feynman gauge is given by —ig^ v D(K) where D(K) is defined in 
(^). The fermion propagator is defined in (0). For the diagram in Fig. [la] we obtain 



where 



G^ a) {P u P 2 ,P 3 ,P,) 



(~,tiv{a) i p p p p \ 
U R2 I 2 ! ' ~2j "31 "ll 



^"("l I p p p p 



^fii/(a) i p p p p \ 
^rara^ (-Pi ) ^2 , ^3 , A ) 

G^W(^1^2,P 3 ,P4) 



(2tt) 4 
rf 4 K 

rf 4 i f 
(2tt, 

rf 4 if 
(27) 

ri 4 iT 
d A K 
d 4 K 



'(if)[aia 2 a 3 a4 + rif 2 f 3 f4 + /ia 2 a 3 a4 
+r 1 f 2 a 3 a i + r 1 r 2 ha i + rif 2 f 3 / 4 ] , 
H^(K)[aia 2 a 3 d4, + rif 2 f 3 f4 + /ir 2 f 3 r 4 
+ai/ 2 a 3 a 4 + axr 2 ha A + air 2 r 3 / 4 ] , 
4 if'"' (if) [01626364 + rif 2 f 3 f 4 + fia 2 f 3 f 4 

+rif 2 r 3 r 4 + aia 2 / 3 a 4 + aia 2 f 3 / 4 ] , 
4 H t * u (K)[a 1 a 2 a 3 a 4 , + rif 2 f 3 f4 + /ia 2 a 3 f 4 
+r 1 f 2 a 3 f i + ri_f 2 ] 3 f^ + aia 2 a 3 / 4 ] , 
if" (if) [(11626364 + rif 2 f 3 f 4 + /ia 2 / 3 a 4 + ri6 4 / 2 / 3 

+rif 2 a 3 a 4 + /ia 2 r 3 / 4 + rif 3 / 2 / 4 + aia 2 f 3 f 4 ] , 
if (.^[(11626364 + rif 2 f 3 f 4 -I- ria 2 a 3 a 4 + aif 2 f 3 fi 
+fif 2 a 3 a 4 + fif 3 f 2 a 4 + A/ 4 r 2 f 3 ] , 



(27T) 



- if" (if) [ai 02030.4 + rir 2 r 3 r4 + a\a 2 a 3 r 4 + r\r 2 r 3 a 4 
+/1/46263 + f 2 f4ria 3 + fahn^] , 



Hfj, u (K) = e A la {If - f i + mb M (f 2 + f 3 + # + ™b^(f 2 + V- + m h* 



(34a) 
(34b) 
(34c) 
(34d) 
(34e) 
(34f) 
(34g) 
(35) 



and we used the notation a\ = Da(K), r\ = Dr(K), fi = Dp(K) for photon propagators and a,j = Da{Kj), 
fj = Dfi(Kj), fj = Dp(Kj) [j = 2,3,4 and K 2 =P 2 +K, K 3 =P 2 +P 3 +K, K 4 =K-P 1 ] for electron propagators. Note 
that the first two terms in each equation of ( |34| ) vanish by contour integration since all poles are located on the same 
side of the real axis. 

Before doing any momentum integrations, we can see immediately from ([34j) that the one- loop retarded functions 
(Cm, i — 1,2,3,4) are linear in the distribution functions, since they are linear in the symmetric propagators fi 
and/or fi. Terms with higher powers of the distribution functions are explicitly cancelled in the Keldysh formalism. 
This type of cancellation also occurs for the 2- and 3-point retarded-advanced functions j27j. In contrast, the mixed 
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retarded-advanced functions are bilinear in the distribution functions. One of the advantages of the Keldysh formalism 
is that it allows us to immediately identify the leading order temperature dependence. 

In the HTL approximation, the external momenta are soft (of order eT with e<l) and the loop momenta are 
hard (of order T) ^,^|. In this approximation, we can neglect the external momenta and the bare electron mass m 
relative to the loop momenta. For H^ U (K) one then obtains 

H^K) « 8e A KnK v # . (36) 

We first consider G"^. Fig. [la] gives 

G^ a \Px,P 2 ,P3,P4) = g J 7^)4 H^(K)[ha 2 d 3 a 4 + rj 2 a 3 a 4 + r 4 f 2 ] 3 a 4 + nr 2 f 3 f 4 ] 

= h+h+h+ h- (37) 



From (|), (@) we have 



Oi(K) =CLi(K) 



K 2 — mf — isgn(fc )e ' 

n{K) = n(K) = — 2 * — , 

K - mf + isgn(fc )e 
fi(K) = -2ni8(K 2 )(l + 2n B (k)) , 

fi{K) = -2iri6(K 2 - m?)(l - 2n F (k)) , (38) 

where i=2, 3, 4 and m 2 —m 3 ^m 4 =m which stands for the electron mass. We calculate the contribution from each of the 
terms in ( |37| ) in the HTL approximation: K ~ T 3> Pi, P%, P 3 , P4,, m. For the temperature-dependent contributions 
we obtain 

h = J dkkn B (k) J dfiVyy • 6336+62 , (39a) 

h = ^3 J dkkn F (k) J dnv^vy -6^636^1 ( 39b ) 
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where 



I 3 = -—— J dkkn F (k) J dClV^VY -63 64623, ( 39c ) 
/ 4 = -t^3 y dkkn F (k) J dQV»V u Y- 6+646+ , (39d) 

6 " = P-T^> (40a) 
^=(^y^' U = 1,2,3,4, (40b) 

and V = (1, V) with V = fe/fc a light-like unit vector, and J dtt denotes integration over the orientation of V. 
By interchanging P 3 and P 4 we obtain the contribution from Fig. [lb] to Gr\: 



Gm b) =Ji + J2 + J 3 + J4, (41) 



with 



Jl = (Sj 1 / dkknB ^ J^v^vy-btb^b^, 



(42a) 



• 4 



J 2 = / dkkn F (k) / dnyyy -63646+2, (42b) 

J3 - 7^ / dk k n F(k) / V»V v y • 1^64 6^ , (42c) 

j 4 = __2JjL / ,//,-/,•„,,(/,••) / r/or''r"v^//; ^7/ t2 . (42<i) 
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Combining (39) and (f42|), the total contribution to G^ x from Fig. [1] reads 

g^(p u p 2 ,p 3 , p 4 ) = GZ {a) + G^ b) = ^ J dn vvy ■ (43) 

where 

e 2 f e 2 T 2 
m 2 =— dkk(n B {k)+n F (k)) = —^- (44) 

is the usual expression for the square of the dynamically generated thermal electron mass. 

Proceeding similarly, we can calculate the other three retarded 4-point functions in finite temperature QED in the 
HTL approximation. The final results, including both diagrams in Fig. [1], are 

ie 2 mi r 

G^ 2 (P 1 ,P 2 ,P 3l P 4 ) = J dQ V^Yb+b+M - bt) , (45a) 

G^n, p 2 , p 3 , p 4 ) = J dn vvYt&b&fc - b-) , (45b) 

G^P^P^P,, P 4 ) = J dQ F7"f6 2 - 3 6 M (^ - b-) . (45c) 

Our results agree with those obtained in the ITF jL8|; to leading order, all of the retarded 4-point functions are 
proportional to T 2 . 

Next we consider the mixed retarded- advanced 4-point functions. We begin with G^ ra . The leading order 
temperature-dependent piece from Fig. [la] is 

G^(P l ,P 2 ,P 3 ,P 4 ) = \ J 0^ H^(K)[ha 2 ha± + ri f 2 f 3 a 4 + / 1 a 2 r 3 / 4 + ^3/2/4] 

= t 1 + t 2 + t 3 + t 4 , (46) 

where the functions Ti ; i = {1,2,3,4} represent the contributions of the four terms respectively. Inserting the 
propagators (pq) and using the HTL approximation we find 



{2^f Cp J ^ V " VU y b t b 2(b 23 - bts) > (47a) 

(^f df) J dnVI±VV y b i^2{bl - H) , (47b) 

(2 4s^ j dQ VVYWtaiK - K) > (47c) 

-dp J dSlV»VYb+b 2 (b+ 2 - b^) , (47d) 



where 



Ti 
T 2 



C/3 = J dkk[n,B(k)(l-nF(k)) - nF(k)(l+riB{k))] , 

dp = -2 / dkkn F (k)(l-n F (k)). (48) 



(2tt) : 



Inserting the standard forms for the Bose-Einstein and Fermi-Dirac distributions we obtain in thermal equilibrium 

c = O, d = -(2hx2)T 2 . (49) 
The contribution from Fig. [lb] to G^ ra is 

= Si + S 2 + S 3 + Si (50) 

with 
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Si 



s 2 = 



(2tt) 3 

e 4 
(2tt) 2 



d/? / dnv»vyb+ 2 b-(b- -b+), 
cp J dnv^yb^b+M -6r), 



5 3 = 



(2tt)3 

,4 



(51a) 
(51b) 
(51c) 
(51d) 



Adding <Wm and (51) we obtain 



' V It )■ 



(27T) 



dfi Wy ^ 6 2 "(&+ (63 + 6 4 ) - 6- 2 (6+ + 6+)) + cp b-(b+(b^ 3 + b 2i ) - + • (52) 



In spite of (MS) we ke ep the terms ~ cp because we will discuss the extension of these results to non-equilibrium 



situations in Sec. VII 



The other two mixed retarded-advanced 4-point functions are similarly calculated: 



(27T) 



(2tt)3 



dn v»v"Y dp b+(bt 2 (bz + hi) - b- 12 (b+ + b+)) + cp 6 2 -(6+(6 2 - 3 + 6 2 - 4 ) - + &2U) 



dn v»v»y cp (b 2 bf(b 23 + 6 2 - 4 ) - b+^(b+ 4 + &+)) 



(53) 



The last function vanishes in equilibrium since cp = 0. The other mixed retarded-advanced 4-point functions are 
again proportional to T 2 . Furthermore, since V 2 =0, all seven 4-point functions are traceless (G£J = 0) in the HTL 
approximation. 



IN THE HTL APPROXIMATION 

-point retarded-advanced 



V. 4-POINT SPECTRAL FUNCTIONS IN HOT QED 

From ([i"3|), (45), (|52] ) and (|53| ) one obtains the following spectral representations for the 4 
functions in the HTL approximation: 

G£j(wi,u;2, u>3, w 4 ) = 7^33 J dQ.idVL-idQ.sdVLi (5(f2i+f2 2 +^3+^4) p^J^i, ^2, ^3, ^4) a^a^a^ - a 2 ) , 

° 4 ) C( fil ' ^2, ^3, ^4) a23 a 24( a r ^ a t) i 



y dOidf2 2 df2 3 df2 4 (5(^1 + ^2 + ^3 + ^4) /9^ L (Oi, 2) ^3, i 
/ d^dttzdttzdni (5(^1 + ^2 + ^3 + ^4) /Cl(°1> fi 2, 03, J 



G^ 2 (wi,w 2 ,W3,w 4 ) = 7-^7- 
( 27I 7 

G^g(wi, W 2 , W 3 , W 4 ) = 1 y d^idri 2 d^3^4 (5(^1+^2 + ^3 + ^4) Phtl(^!' ^3, ^4) a23 a 24( a r _ a 2 



(54a) 
(54b) 
(54c) 
(54d) 



G^(U)1, UJ 2 , U} 3 , W4) = T^Tg yd^ldfi2^3^4 (5(^1+^2 + ^3 + ^4) Phtl(^1' ^2, ^3,^4) a 2 3 a 24( a l 

G{£ ro (wi, «3, w 8 , w*) = 7A3 y dfl^dQadQi <5(n 1 +0 2 +n 3 +0 4 ) p£ L (fii, n a , fi 3 , fi 4 ) 

ai a^a^ag + a 4 ) - a^a 4- + a t )) + "2 (% ( a 23 + fl 24) - a 2~ a r( a 24 + a 23)) 

G^, r (wi,wa,w3,w 4 ) = 77^3 y dn^dQadQi < s(n 1 +n 2 +n 3 +n 4 ) ^ L (ni,n 2 , n 3 , n 4 ) 

a l"( a 23 + °24) - a 2~ a r( a 24 + ^3) ) ( 54r ) 



(54e) 



(012(03" + 4 ) - or2( a 3 + O4 )) + a 2 a 2~ a l 

G^jwi,^,^,^) = 77^3 y dn 1 dn 2 dQ 3 dn i 8(sii+si 2 +sh+^) pZl (°i> ^2,n 3) o 4 ) 

■ "2 (a2~ a l~( a 23 + «24) - a 2 a l( a 24 + a 23)) i 

where 



(54g 
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p^ TL (n 1 ,n 2l n 3l n 4 ) = 2Tr 2 e 2 mf j f dnv»v v Y5(P 1 -v)5(p 2 -v)5(P 3 -v) (55) 

(with Pf = tii, Pi = Pf, i={l, 2, 3} on the r.h.s.) and 

dpe 4 cpe 4 
Ql = o- 2 2 2 i "2 = . 2 2 2 ■ (56) 



These equations also hold in the non-equilibium situations studied in Sec. VII. In equilibrium we use (49) to obtain 



■ 81n2 n fx* 

ai = i — — ; a 2 = (57) 

Note that the frequencies fi, should not be confused with the angular factor fl. From V 2 =0 it follows that 

OW)£ = 0. (58) 



In contrast to the general spectral representations (gOj) of the 4-point functions given in Sec. |HI], in the HTL 
approximation the spectral representations ( |54| ) involve only a single, real- valued spectral density P^ L - This agrees 
with the result of Taylor Jl8| who showed for QCD in the ITF that in the HTL approximation all spectral densities 
degenerate to a single function. Clearly the spectral density (p5|) is symmetric in the Lorentz indices and under 
interchange of the external momenta. Using the 5-functions in (|55[) it is also easy to show that is transverse with 
respect to all external momenta: 

p >ht L =°> i = l,2,3,4. (59) 
A previous RTF calculation of the 3-point spectral function in hot QED in the HTL approximation gave J2(| 

iremi f 

^ TL (n 1 ,fi 2 ,fi 3 ) = — / dQV^Y 6{Pi-V)5(P 2 -V). (60) 
By comparing (|55|) with ( |59| ) we find the "sum rule" 

dn 3 P f TL (n u n 2 ,n 3 , -(fi 1+ n 2 +n 3 )) = 4^ e ^ TL (o x , n a , -(n x +n 2 )) . (6i) 

This agrees with the ITF results in E^j. A similar sum rule exists between the 2- and 3-point spectral densities: the 
2-point spectral density is given by 

,w(fii,-ni) = ^ / dny*(Pi-vo, (62) 

and the temporal component of the 3-point spectral density /?h TL (^i> ^2, — (^1+^2)) obeys the sum rule 

J d0 2 p° TL (f2 1) Q 2) -(n 1 +0 2 ))=7re / 9 HTL (fi 1) -fi 1 ). (63) 
We believe that analogous sum rules exist for the higher order n-point functions. 

VI. WARD IDENTITIES BETWEEN 3- AND 4-POINT FUNCTIONS IN THE HTL APPROXIMATION 

In this section, we verify the HTL Ward identities between the 3- and 4-point vertex functions in QED, without 
doing momentum integrations. We split propagators in the integral expressions (|34|) by using the identities p8| 

Dr/a(K 2 ) D r/a (K 3 ) = 1 (Dr/a{K 2 ) - D R/A (K 3 )) (64) 
This trick allows us to rewrite (|34|) as 
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^Rl 



Kjr R2 



R3 
T R4 

vara 

r<fJ,v(a) 
rraa 

- rn.n.r 



d 4 K AK^K"^ 
(2tt) 4 2P 3 -K 
d 4 K <±K^K v Ij: 
(2tt) 4 2P 3 -K 
d 4 K <±K^K v Ij: 
(2tt) 4 2P 3 ■ K 
d 4 K 4K^K" fC 
(2tt) 4 2P 3 -K 
d 4 K AK^K^If 
(2tt) 4 2P 3 ■ K ^ 
d 4 K AK^K V I£ 



[/i(fl2 - 03)04 + rif 2 a4 - r 4 f 3 d 4 + ri(f 2 - f 3 )/ 4 ] , 
[fi(r 2 - r 3 )f 4 + aif 2 d 4 - aif 3 a 4 + ai(f 2 - f 3 )/ 4 ] , 
[fi(a 2 - f 3 )r 4 + ri/ 2 f 4 - ai/ 3 a 4 + ai(a 2 - f 3 )f 4 ] , 
[/i(«2 - 03)^4 + rif 2 f 4 - r 1 f 3 f 4 + oi(a 2 - a 3 )/ 4 ] , 
-fihdi + r 1 (f 2 - d 3 )d 4 + /i(a 2 - f 3 )/ 4 + ri/2/4 + ai(a 2 - f 3 )f 4 ] , 



(2tt) 4 2P 3 ■ K 
d 4 K AK^K v If 
(2tt) 4 2P 3 -# 



[ri(a 2 - a 3 )a 4 + ai(r 2 - r 3 )r 4 + /i/ 2 a 4 - /i/ 3 a 4 + fif 4 (r 2 - r 3 )] 
[oi(a 2 - a 3 )f 4 + n(f 2 - f 3 )a 4 + A/ 4 (a 2 - a 3 ) + / 2 / 4 ri - / 3 / 4 ri] 



where we have used the HTL approximation to write K 3 — K\ = 2K ■ P 3 + P 3 « 2if • P 3 



P1+P2 




K-Pl 



P2 



(65a) 
(65b) 
(65c) 
(65d) 
(65e) 
(65f) 
(65g) 



Fig. [2]: 3-point vertex in QED 



The integral expressions for the 3-point vertex functions in QED are obtained from Fig. [2]. In the Feynman gauge 
we have 



G R (P 1 ,P 2 ,P 3 +P 4 ) 


= e 3 


G Ri (Pi,P 2 ,P3- 


KP 4 ) 


= e 3 


G Ro {P u P 2 ,P 3 - 


hP 4 ) 


= e 3 


G»(Pi,P 2 ,P 3 - 


hP 4 ) 


= e 3 


Gpi(Pi, P 2 , P3- 


hP 4 ) 


= e 3 


G» (P U P2,P 3 - 


hP 4 ) 


= e 3 



(2^) 4 
ri 4 i f 
(2tt 
d 4 K 



{2K >1 I^)[aia 2 a 4 +r 1 r 2 r 4 +f 1 r 2 r 4 +aif 2 a 4 +a 1 r 2 f 4 \, 
4 (2KVtf)[aid 2 a 4 +r 1 f 2 f 4 +f 1 a 2 d 4 +r 1 f 2 a 4 +r 1 r 2 f 4 ], 



(2n) 4 
d 4 K 

W) 
d 4 K 



(2K >1 I^)[aia 2 a 4 +r 1 r 2 r 4 +f 1 a 2 r 4 +rif 2 r 4 +aia 2 f 4 \, 
4 {2K^If)[aid 2 f 4 +rif 2 d 4 +rif 2 f 4 +fia 2 f 4 }, 



(2ir) 4 
d 4 K 
W) 



(2K tl I^)[a 1 r 2 a 4 -\-r 1 a 2 r 4 +aif 2 f 4 +fif 2 r 4 }, 
4 (2K^I^)[ria 2 a 4 +aif 2 f 4 +f 1 f 2 d 4 +fif 2 f 4 \, 



(66a) 
(66b) 
(66c) 
(66d) 
(66e) 
(66f) 



where for later convenience we have chosen to write on the l.h.s. — (P1+P2) = P3+-P4. Comparing ( |65[ ) and ( |6q ) we 
obtain the following Ward identities between the vertices shown in Fig. [la] and Fig. [2]: 
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P^G^ a) (P l .P 2 , P 3 , P 4 ) = e(G v m {P 1 ,P 2 , P3+P4) - G" m (P 1 ,P 2 +P 3l P i )) , (67a) 

P^GZ ia) (Pi, P2, Ps, Pi) = e(G' J R (P 1 ,P 2 ,P 3 +P 4 ) - G R (P U P 2 +P 3 , Pa)) , (67b) 

P^G^ a \P 1 ,P 2 ,P 3 , Pi) = e(G v Ro {P x ,P 2 ,P 3 +Pi) - G v R {PuPi+P3,Pij) , (67c) 

P^G$ a \P u P 2 ,P 3 , Pi) = e(G Ro (P u P 2 ,P 3 +P i ) - G^jPi.P^P^Pi)) , (67d) 

^G^)(P 1 ,P 2 ,P 3! Pi) = e(G» F (P u P 2 ,P 3 +P 4 ) - G Fo (P u P 2 +P 3 , P 4 )) , (67e) 

P 3m G^(Pi,P 2 ,P 3 , Pi) = e(G F (P ll P 2l P 3 +Pi) - GpiPuP^P^Pi)) , (67f) 

Ps^G^OPi'^Ps, P 4 ) = e(G Fo (P u P 2 , P 3 +Pi) - G Fo (P u P 2 +P 3 ,Pi)) . (67g) 
Similarly one can obtain the Ward identities satisfied by the diagrams in Fig. [lb] and Fig. [2]: 

P 3lJ ,G R f\P x ,P 2 ,P 3 ,Pi) = e(G Rt (P 1 +P 3l P 2l P 4 ) - G Rl {P u P 2 ,P 3 +Pi)) , (68a) 

P3»G% 2 {b) (P 1 ,P 2 ,P 3 ,Pi) = e(G R {Pi+P 3 ,P 2l Pi) - G" R (Pi, P2, P3+P4)) , (68b) 

P^G^ b) {P u P 2 ,P 3 ,Pi) = efa^+P^P^ P 4 ) - G Ro {P u P 2l P 3 +Pi)) , (68c) 

P^G^ h) (Pi,P 2 ,P 3 ,Pi) = e(G Rl (P 1 +P 3 ,P 2 ,Pi) - G Ro {Pi,P 2 ,P 3 +Pi)) , (68d) 

P-^GZi h XPu Pi, Ps, Pi) = e(G^(P!+P 3 , P 2 , P 4 ) - G F (P U P 2 ,P 3 +Pi)) , (68e) 

P^G^) [Px,P 2 ,P 3 ,Pa) = -eG F {P 1: P 2: P 3 +Pi) , (68f) 

P^G^iPuP^P^Pi) = e(G Fo (P 1 +P 3 ,P 2 ,Pi) - G Fo (P u P 2 ,P 3 +Pi)) . (68g) 

By combining ( |67| ) and ([38]) one obtains the following HTL Ward identities between the 3- and 4-point vertex functions 
in QED in the RTF: 

P^G^{P U P 2 , P 3 , P 4 ) = e(G» m (Pi+P3,P2, Pi) ~ G Rz (P l7 P 2 +P 3l Pi)) , (69a) 

P^G^Px, P 2 , P 3 , P 4 ) = e(G R {P 1 +P 3 ,P 2 ,P i ) - G R {P u P 2 +P 3 .Pi)) , (69b) 

P^G^Pu P 2 , P 3 , P A ) = e(G R (P 1 +P 3 ,P 2 ,P 4 ) - G v Ro (P u P 2 +P 3 ,P 4 )) , (69c) 

P 3t iG R u 4 (P ll P 2 ,P 3 ,Pi) = e(G Ro {Pi+P 3l P 2 ,Pi) - G Ri (P u P2+P3, Pa)) , (69d) 

P^G% ra {Pi, P 2 , P 3 ,P 4 ) = e^GpiPi+P^P^Pi) - Gp^Pi.P^P^Pi)) , (69e) 

P^GZaA p u P2, P 3 , Pa) = -eG F (P u P 2 +P 3 , Pi) , (69f) 

P^G^Px, P 2 , P 3 , Pi) = e(G Fo (P 1 +P 3 ,P 2 , Pi) - Gp^PuP.+P^Pi)) . (69g) 

Notice that in equilibrium the last equation gives simply = (see (E9I) and (53)). These results are structurally 



similar to the zero temperature Ward identity |29J and to the Ward identities between the retarded 3- and 4-point 
HTL vertices in QCD obtained within the imaginary time formalism ||] or from kinetic theory Jl0| . 

VII. EXTENSION TO NON-EQUILIBRIUM SITUATIONS 

The HTL method has been widely used in qualitative studies of the QGP and in other explicit calculations at finite 
temperature [^). This powerful method was derived within the ITF for equilibrium field theory, and this is where its 
range of applicability can be clearly defined |4|,||,|| : HTL resummation is required and applicable for the study of soft 
processes in a weakly interacting plasma, with momentum scale p^gT^T(g^l) much below the typical thermal 
or "hard" momenta in the medium. Realistic physical systems, on the other hand, are frequently out of equilibrium, 
which means that calculations must be carried out in the RTF. As long as the local momentum distribution f(x,p) 
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doesn't deviate too strongly from a thermal one, and can still be roughly characterized by a parameter T which gives 
the scale for the typical momenta at point x, a similar "hard loop" (HL) resummation scheme should be applicable for 
the study of the dynamics of soft modes (p ~ g T <C T) in a weakly interacting (g <C 1) non-equilibrium system ps[ . 
This approach is expected to be reliable for equilibration processes requiring momentum exchanges which happen 
slowly, on a time scale r ~ (g T) ^> 1/p . This reasoning does not apply to arbitrary non-equilibrium situations, 
but only to plasma states which are sufficiently close to global thermal equilibrium; on the other hand, deviations 
from chemical equilibrium can be arbitrarily strong without invalidating the scheme. 

In this section we discuss the generalization of our results to this particular type of non-equilibrium scenarios. 
Since we have used the Keldysh or closed time path contour, the procedure is straightforward. Equations (|l|) and 
(||) can be used in non-equilibrium situations by replacing the x-independent equilibrium Bose-Einstein and Fermi- 
Dirac distributions by non-equilibrium distributions (Wigner functions) fs(x,p) and fp{x,p) which depend on the 
space-time coordinate and the four momenta (21[ . Equations (^) and (|t]) remain valid for bare propagators, with the 
equilibrium distribution functions (nB,np) replaced by the non-equilibrium ones (/b,/f), but the last equations in 
each of (||) and (|), (i.e. the fluctuation-dissipation theorem which reflects the KMS condition), no longer hold for 
full propagators | 2q| . 

Equations ( |30| ) give the spectral representations of the seven 4-point functions ( ]io| ) in equilibrium. These expressions 
remain valid out of equilibrium since the KMS condition was not used in their derivation. Out of equilibrium there 
are eight more 4-point vertex functions (G ar rr, G ra rr, G rrar , G rr ra, G ar ar, G aa rr, G ar ra and G rr rr) with similar 



spectral representations. They can be derived using the same methods as those given in Sec. [II. 

The results in Sees. IV and [v] can be generalized to non-equilibrium situations by simply replacing in ( |44| ) and ( |4^ ) 
the Bose-Einstein and Fermi-Dirac distributions w_b(p°) and np(p°) by non-equilibrium Wigner functions Jb{ x tP) an d 
fp(x,p). The thermal mass mp in ( [44] ) and the factors cp, dp in (^) become functions of x, m(x), c(x), d(x), and it is 
only in the special case of thermal equilibrium that c = as obtained in ( |49| ) . These generalized results correspond to 
the "hard loop" (HL) approximation in slightly off-equilibrium plasmas and are applicable in the situations discussed 
at the beginning of this section. We have shown that the 7 components of the 4-point function given in ( |To| ) correspond 
to 5 spectral densities which degenerate in the HL approximation to a single real function. This result remains valid 
out of equilibrium, where the KMS conditions do not hold. The eight additional 4-point functions which must be 
included out of equilibrium will involve more spectral functions, but in the HL approximation it is straightforward 
to show that they degenerate to the same HTL spectral density (p5|), with the thermal mass mp replaced by the 
generalized mass m(x). 



The one-loop Ward identities in Sec. VI were obtained without explicitly doing any momentum integrals; the result 



thus holds independently of the form of the distribution functions (which enter through the functions fi, fi in (65) and 
(|66|)). Out of equilibrium there will be additional Ward identities for the eight additional vertex functions which can be 
derived in the same way. For example, the Ward identities that correspond to the three functions G aa rr , G ar ra , G a 
have been calculated, and they have the same form as those obtained in (69 e,f,g). 



T ara r 



VIII. CONCLUSIONS 

We have studied the 4-point vertex functions in the Keldysh formulation of real-time thermal field theory |^2| . This 
formalism has recently gained increased popularity because it avoids the need for analytical continuation associated 
with the imaginary time formalism, and it allows for a generalization to non-equilibrium situations. Starting from the 
largest time and smallest time equations which relate the 16 components of the 4-point real-time vertex function, we 
have derived spectral integral representations for the 7 retarded-advanced 4-point functions ([l0|). We have explicitly 
calculated these 7 functions for hot QED at 1-loop order in the HTL approximation, using the Keldysh representation. 
In particular, we have demonstrated the usefulness of this representation for the calculation of many-point functions 
and showed how terms with high powers of the distribution functions cancel explicitly before any momentum inte- 
grations are done. The fully retarded vertex functions are found to be linear in the distribution functions p7j ]. In the 
HTL approximation, they have the same structure as their QCD counterparts which were derived in the ITF |^,^0| as 
well as from kinetic theory [|l0| . We have also calculated the three mixed retarded-advanced 4-point functions which 
are needed for a complete description of the real-time 4-point vertex and which have not previously appeared in the 
literature. They are quadratic in the distribution functions. One of them (G rraa ) was found to vanish in thermal 
equilibrium, while the other 4-point functions are all proportional to T 2 . 



We have derived spectral integral representations for the 7 retarded-advanced 4-point functions (10). The spectral 
densities were calculated explicitly for hot QED in the HTL approximation and found in this case to degenerate to 
a single real spectral function. Its trace over the Lorentz indices vanishes, and it is transverse with respect to all 4 
external momenta. It obeys a simple sum rule with the correspondingly approximated 3-point spectral density. 



14 



By contracting the 7 retarded-advanced 4-point functions (|10|) for QED in the HTL approximation with one of the 
photon momenta we have derived one-loop Ward identities. This calculation was done without doing any momentum 
integrations, by comparing the integrands of the contracted HTL 4-point vertices with those of the HTL 3-point 
vertices. The resulting set of real-time Ward identities at finite temperature generalizes the zero temperature Ward 
identity and can be compared with the HTL Ward identity from the ITF and from kinetic theory. Due to the matrix 
form of the real-time thermal Green functions, the Ward identities have a more complex structure: one finds a whole 
class of finite temperature Ward identities which relate retarded-advanced 4-point functions to retarded-advanced 3- 
point vertex functions. If either the ingoing or the outgoing fermion leg has the largest time, the Ward identity involves 
only one kind of retarded-advanced 3-point vertex functions, and the right hand side vanishes when the contracted 
photon momentum (P3) goes to zero. If one of the photon legs has the largest time, two different retarded-advanced 
3-point functions are involved, and the right hand side no longer vanishes (Eqs. (|69|c-f); the corresponding component 
of the 4-point vertex is thus singular in the zero-momentum limit P3 — > |pl| . Similar features were found for the 
Ward identities relating the HTL 3- and 2-point functions in hot QED |^0| . 

Our Ward identities are not general, in the sense that they were not directly derived by functional methods from 
the RTF path integral. However, in the case of the 3-point functions it is known that the one-HTL Ward identities 
p0[ are identical in structure with those satisfied by the exact vertex functions |52) . We therefore expect the same to 
be true in the present case. 

Finally, we have discussed the generalization of our results to certain classes of "moderately non-equilibrium" 
situations (Sec. VII) in which a generalized HTL approximation ("hard loop approximation") can be applied. Out of 
equilibrium the number of independent 4-point functions is larger, but in the HL approximation the spectral functions 
have the same structure as in equilibrium. All of the nice features of equilibrium HTLs persist in such "moderately 
non-equilibrium" situations. 
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I. INTRODUCTION 

One of the most significant advances in quantum field theory at finite temperature [1-3] is the effective perturbation 
theory of Braaten and Pisarski [4,5] which is based on the resummation of hard thermal loops (HTLs). HTLs were 
originally obtained by computing one-loop diagrams in the imaginary time formulation of finite temperature field 
theory. HTLs are ultraviolet finite, gauge invariant and satisfy simple Ward identities [6-9]. These remarkable 
properties have been traced to the fact that HTLs describe simple semiclassical physics [6,10,11]. In particular, the 
HTLs can be derived from a set of classical kinetic equations [6,10] as well as from a non-local effective Lagrangian 
[8]. The computation of HTLs is fairly technical because of their non-trivial momentum and energy dependence, but 
can be simplified by using the Ward identities [4]. 

The analytical structure of Green functions at zero or non-zero temperature, in the imaginary or real-time formula- 
tions, is controlled by their spectral densities, combined with asymptotic boundary conditions. At finite temperature 
the spectral functions can also be directly related to transport coefficients [12,13]. The study of spectral functions 
helps us to understand the quasi-particle structure of a field theory and to identify the microscopic processes underly- 
ing the dynamics. But Green functions and their spectral functions are not easy quantities to evaluate perturbatively 
at non-zero temperature, especially for many-point functions. The two-point functions and their spectral densities 
have been widely studied and applied to quark gluon plasma (QGP) investigations [14,3,5,12]. Our knowledge about 
many-point spectral densities is still much less comprehensive [15-20]. The spectral representations of 3-point func- 
tions for self-interacting scalar fields were discussed in Refs. [15,16,19,20]. In [18] 3- and 4-point spectral densities for 
pure gluon dynamics were calculated in the HTL approximation using the imaginary time formalism (ITF) [4], and 
in the same approximation the 3-point spectral density in QED has been calculated within the real time formalism 
(RTF) [20]. 

The popularity of the RTF is diminished by the technical complications resulting from the doubling of degrees of 
freedom [21]. However, in the ITF, one has to perform an analytic continuation of the imaginary external energy 
variables to the real axis, which is avoided in the RTF. In addition, the ITF is restricted to equilibrium situations, 
while the RTF can be extended to investigate non-equilibrium systems [21-23]. 

In this paper we study 4-point functions in QED using the real-time formalism. We adopt the Kcldysh or closed 
time path (CTP) contour [22]. In Sec. II we introduce the Keldysh representation. The spectral representations of 
the 4-point functions is derived in Sec. III. In Sec. IV we use the Keldysh representation to calculate the 4-point HTL 
vertex functions in hot QED. The 4-point spectral densities in the HTL approximation are extracted in Sec. V and 
shown to degenerate to a single real spectral function. In Sec. VI we derive real-time Ward identities between the 3- 
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and 4-point HTLs and their spectral functions. The generalization of our results to certain types of non-equilibrium 
situations is discussed in Sec. VII, while Sec. VIII gives a summary of our results. 



II. KELDYSH REPRESENTATION 

We will use the Keldysh representation [21,22] of the real-time formalism (RTF) throughout this paper. In this 
representation the single-particle propagator for free bosons has in momentum space the form [21,22] 



D(K) = [ **=^ \ 27ri(5( ^ 2 _ m2) ( Mfco) 0(-k ) + n B (k ) j (1) 

V K*-m*-ie J J \ 6(k ) + n B (k ) n B (k ) J W 

where K = (k ,k), k = \k\, 9 denotes the step function, and the equilibrium distribution function is given by 
n B (k ) = l/[cxp(|fc |/T) — 1]. For fermions the bare propagator can be written as 



S(K) = (% + m) 



g 2_ m 2„ le / ' \9(k ) -n F (k ) n F (k ) 



(2) 



where the Fermi-Dirac distribution is given by n_p(fco) = l/[exp(|fc |/T) + 1]. The components of these propagators 
are not independent, but fulfill the relation 

Gn — G12 — Gii + G22 = 0, (3) 

where G stands for D or S, respectively. 

By an orthogonal transformation of these 2x2 matrices one arrives at the representation of the propagators in 
terms of advanced and retarded propagators which was first introduced by Keldysh [22]. The three nonvanishing 
components of this representation are [21] 

Gr = Gn — Gu, 
Ga = Gn — G21, 

G F = Gu+G 2 2. (4) 



The inverted relations read 



G11 = 2 (Gf + Ga + Gr), 
G12 = — (G F + Ga — Gr), 
G21 = 2 {Gf — Ga + Gr), 

G22 = 2 (Gf - Ga - Gr). (5) 



The relations (3-5) also hold for full propagators. 

Using (1) and (2) in (4), the free propagators are given in the Keldysh representation by 

Dr(K) : 



K 2 — m 2 + i sgn(/c )e 

D A (K) 



K 2 — m 2 — zsgn(/c )e ' 

D F (K) = (l+2n B (k )) sgn(fc ) (d r (K) - D A (K)^j = -2m (l+2n B (k )) S(K 2 -m 2 ) (6) 



for bosons and 



- -77-2 H m riH = + m)D R (K) . 

K 2 — m z + i sgn(fc )e 

S a(K) = -77-2 f + m 77-7 =(% + m)D A (K) . 

K 2 — m 2 — i sgn(fc )e 



S F (K) = (l-2n F (fc )) sgn(fco) (s R (K) - S A (K)) 

= -2m (If + m) (l-2n F (fc )) S(K 2 - m 2 ) = (If + m)D F (K) (7) 
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for fermions, where the last equation in each of (6) and (7) is a consequence of the dissipation-fluctuation theorem. 
Although Da.r = Da,r we keep the bar in our notation to more easily recognize fermion and boson propagators in 
the calculations below. 



III. SPECTRAL REPRESENTATION OF THE 4-POINT VERTEX 



In this section we briefly review some useful relations among the different thermal components of the 4-point 
functions, and we derive their spectral representations. Some of these relations have already been reported in the 
literature [14-16,24] using different notation. We discuss these results at the end of this section. For any 1PI four- 
point function we define the four incoming external momenta as Pi, P2, P3, and P 4 = — (P1+P2+P3). In this section 
we suppress all indices except Keldysh indices. 



A. Largest and smallest time equations 

If £4 is the largest time argument, one can obtain the "largest time equations" [14,15,24]: 

#43 #32 #21 (Gabcl + G a bc2) = #43 #31 #12 {Gabcl + G abcl) = , (8a 

#42 #23 #31 {Gabcl + G a bc2) = #42 #23 #13 {G a bcl + G a bc2) = , (8b 

#41 #13 #32 {G a bcl + G a bc2) = #41 #12 #23 {G a bcl + G a bc2) = , (8c 

#32 #21 #14 {Gable + G a b2c) = #31 #12 #24 {G able + G a b2c) = , (8d 

#23 #34 #41 {Galbc + G a 2bc) = #24 #43 #31 {G a lbc + G a 2bc) — , (8e 

#13 #32 #24 {Glabc + G2abc) = #12 #23 #34 {Glabc + G2abc) = , (8f 

where a, 6, c can be either 1 or 2, and Oij = 9{U — tj) is the step function. 

By tilde conjugation (which replaces time ordering with anti-time ordering [21,23]) one obtains equations of the 
form: 

#12 #23 #34 {Gabcl + G a bc2) = #13 #32 #24 {G a bcl + G a bc2) = , (9a) 

#41 #12 #23 {Gable + G a b2e) = #42 #21 #13 {G a blc + G a b2c) = , (9b) 

#14 #43 #32 {Galbc + G a 2bc) = #13 #34 #42 {G a lbc + G a 2bc) = , (9c) 

#42 #23 #31 {Glabc + G2abc) = #43 #32 #21 {Glabc + G2abc) = , (9d) 

Eqs. (8) and (9) are the analogues of the "largest time equations" and "smallest time equations" , respectively, of 

Ref. [14]. They will be used extensively in the derivation of the spectral representations below. Their generalization 
to arbitrary n-point functions is straightforward. 



B. Derivation of spectral representations 

One can construct the "retarded-advanced" vertex functions from the sixteen components of the real-time 4-point 
function as in [19,21]. The KMS conditions allows to reduce them to the following seven combinations [25,26]: 

(10a) 
(10b) 
(10c) 
(lOd) 
(10c) 
(lOf) 
(10g) 

It is straightforward to show that the first four are the standard retarded vertex functions 

Graaa — G Rl J G a raa ~ G R2 j G a ara = G R3 J G aaar = G m . (11) 



Graaa 


= Gini 


+ G1211 


+ G1121 


+ G1112 


+ G1122 


+ G1221 


+ G1212 


+ G1222 


Garaa 


= Gini 


+ G2111 


+ G1121 


+ G1112 


+ G1122 


+ G2121 


+ G2112 


+ G2122 


Gaara 


= Gim 


+ G2111 


+ G1211 


+ G1112 


+ G4212 


+ G2211 


+ G2112 


+ G2212 


Gaaar 


= Gini 


+ G2111 


+ G1211 


+ G1121 


+ G1221 


+ G2211 


+ G2121 


+ G2221 


Grara 


= Gini 


+ G1211 


+ G1112 


+ G1212 


+ G2121 


+ G2122 


+ G2221 


+ G2222 


Grraa 


= Gini 


+ G1112 


+ G1121 


+ G1122 


+ G2211 


+ G2212 


+ G2221 


+ G2222 


Graar 


= Gini 


+ G1121 


+ G1211 


+ G1221 


+ G2112 


+ G2122 


+ G2212 


+ G2222 
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The last three we will call "mixed retarded-advanced" functions. 

There are eight other vertices which one can obtain from (10) using the KMS conditions [25,26]. We write down 
only the three which we need: 

Garar = Gnu + Gn 2 i + G2111 + G2121 + G1212 + G1222 + G2212 + G2222, (12a) 
Garra = Gnu + Gm2 + G2111 + G2112 + G1221 + G1222 + G2221 + G2222, (12b) 

Gaarr = Gnu + Gl211 + G2IH + G22II + Gn22 + G1222 + G2122 + G2222- (12c) 

We can re-express these vertices in terms of spectral functions by making use of various properties of theta functions. 
We will use 

#12 + #21 = lj 012021 = 0, #23#34#24 = 023^34- (13) 

(i) We begin with Gri. We rewrite Gri as 

Gri — (023 + #32)(#34 + #43)(#24 + ^42)G_R1. (14) 



Inserting the identities 



into (14) we obtain 



Using the identities 



h^Gm — OuGri — OuGri — Gri (15) 



Gr\ — 012#23#34Gi{i + 0i2024043G_Ri + 614642 623 Gri 

+ 6i404 3 9 S 2Gri + #13#34#42G_R1 + 013#32#24 Gri . (16) 



9(iijk)G R2 = 9 li e ij e jk G R2 = 0, 

9(lijk)G m = 9 u 9 ij 9 jk G RS = 0, 

9(lijk)G R4 = 9 li 9 ij 9 jk G Ri = 0, (i,j,k = 2,3,4) (17) 

which result from conflicting ^-functions, we have 

G R1 = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432)] (G m - G R2 + G R3 - G BA ) 

= [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432)] p 12 (18) 

with 

Pl2 = Gri — Gr2 + Grz — Gr 4 (19) 
= Gi211 — G2122 + 2(Gi212 — G2121) + G1112 — G2221 + G1222 — G2111 + G2212 — G1121. 

(ii) For Gr2, Gr$ and Gr4 one proceeds similarly and obtains 

G R 2 = [0(2134) + 0(2143) + 0(2314) + 0(2341) + 0(2413) + 0(2431)](-pi 2 ), (20a) 
G R3 = [0(3124) + 0(3142) + 0(3214) + 0(3241) + 0(3412) + 0(3421)] p 12 , (20b) 
G R 4 = [0(4123) + 0(4132) + 0(4213) + 0(4231) + 0(4312) + 0(4321)](-pi 2 ). (20c) 

(iii) The procedure for the mixed retarded-advanced four-point functions is similar. We write 

Grraa = (012 + 2 l)(014 + 04l)(023 + 032)(034 + 043)(#24 + 042)G rraa (21) 

and use 

9(3ijk)G rraa = 0, i,j,k = 1,2,4, 

9{Aijk)G rraa = 0, i, j, k = 1, 3, 4, (22) 
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to obtain 

Grraa = [6»(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(2134) + 0(2143) + 0(2314) + 0(2341) + 0(2413) + 0(2431)] G rrac 

Using the identities 

9(lijk)G rraa = 0, i,j,k = 2,3,4, 
6{2ijk)G rraa = 0, = 1,3,4 

allows us to write 

G rraa = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(2134) + 0(2143) + 0(2314) + 0(2341) + 0(2413) + 0(2431)] p 3 

with 

P3 G rraa G rraa . 

Analogously, we can derive 

G rara = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(3124) + 0(3142) + 0(3214) + 0(3241) + 0(3412) + 0(3421)] p 4 

G raar = [0(1234) + 0(1243) + 0(1324) + 0(1342) + 0(1423) + 0(1432) 

+0(4123) + 0(4132) + 0(4213) + 0(4231) + 0(4321) + 0(4312)] p 5 

where 

Pa G rara G rara 
Pb G raar G raar 

Using the Fourier integral representation of the function, 

1 r°° p-in(u-tj) 



13 ' 2niJ_ 00 Cl + ie ' 
it is straightforward to derive the spectral integral representations in momentum space: 

Gfll(wi,W2,W3 1 W4) = j> at[ai2( a 3 + a 4 ) + a 13( a 2 + a 4 ) + a li( a 2 + a 3~)] Pl2 i 

G R2 (uji,uj2,u 3 ,uj4) = j> a\ [a2i( a 3 + a i) + a 23( a r + a i) + <4i( a r + a 3 )]{-Pl2) 
G m {wi, uj 2 , w 3 , w 4 ) = <j> a^a^a^ + a 4 ) + a% 2 {a{ + aj)a^" 4 (aj" + )] P12 , 
Gfl 4 (wi, ^2,^3,^4) = j> a^aj^a^ + ajf ) + a^ 2 (a^ + ) + a| 3 (af + a 2 ")](-pi 2 ) 

G rraa (wi, ^2,^3,^4) = <j> (at[ a 12( a 3 + a 4 ) + a 13( a 2 + a 4 ) + a 14( a 2 + a 3~)l 

+ a 2 [ a 2l( a 3 + a 4 ) + a 23( a r + a 4 ) + a 24( a r + a 3~)l) P3 , 

G™ra (wi, w 2 , u 3 , w 4 ) = ^ (a^a^dg + aj) + a^a^ + aj) + a^a^ + O3 )] 

+«3 [ a 3l( a 2 + a 4 ) + 032 (°r + a 4 ) + a 34( a r + a 2~)]) Pi 
G raor (wi,W 2 ,W 3 ,u; 4 ) = ^ (ai"[ai" 2 ( a 3 + a 4 ) + fl 13( a 2 + a 4 ) + a 14( a 2 + a 3~)l 

+a+[al l (a 2 + ) + a| 2 (af + ) + aj 3 (a^[ +a 2 )]j p 5 



where 



1 



fli ± ie 
1 



flj ± ie 



i,j = 1,2,3,4. 



(31a) 
(31b) 
(31c) 



The frequency arguments of the spectral functions under the integrals are ^(Sli, Sl 2 , ^3, O4) with f2i+f2 2 +f2 3 +f24=0. 
The spatial momenta P\,p 2 ,Pz, and P4—~ (P1+P2+P3) are the same on both sides of these equations and have 
therefore been suppressed. 

We have written the seven vertex functions (10) in terms of one complex spectral density (pi2 = p\ + ip-i) and 
three purely imaginary spectral densities (^3, P4, P5). This means that all sixteen 4-point functions can be expressed 
in terms of five real spectral densities: 



pi = Re [n B 1 (p 20 )G 2122 + 2n B 1 (p 20 + P4o)G 2 i 2 i + n B 1 (p 40 )G 2 22i - n B 1 (p 10 )G 12 22 - n^ip^G 2212} , 
p 2 = -Im[n F 1 (p2o)G 2 i22 + 2n F 1 (p 2 o +P4o)G 2 i2i + n F 1 (p m )G 2 22i - n F 1 (p w )G 12 22 - n F 1 (p 30 )G 2212] , 

P3 = -iP3 = 2 Im G r raa , 
Pi = -iPi = 2 Im G rara , 
P5 = ~iP5 — 2 Im G r aar , 

where we used the following relations in momentum space resulting from the KMS condition [21]: 



G 12 ii = e* 20 G^ 122 , 
G2111 = e^ 10 G^ 222 , 

Grraa — ^ rraa ' G T ara 



G1212 
G1121 

if 



2121 ' 



„/3(P20+P40) (~<* 

e Lz 2212 i 

roar — , 



G1112 = e^ 40 Gt 



2221 ' 



(32a) 
(32b) 
(32c) 
(32d) 
(32e) 



(33a) 
(33b) 
(33c) 



These results should be compared with the expressions derived in [15] in which the spectral representations for the 
four retarded functions Gm, i = 1, 2, 3,4 were given in terms of three real spectral densities. A relation between these 
three spectral densities exists which can be employed to reduce the number of independent spectral functions to two; 
this is consistent with our result. Spectral representations for the mixed retarded-advanced 4-point functions have 
also been previously discussed: in [17], using a different definition of the mixed retarded-advanced 4-point functions, 
the spectral representation for one of these vertices was given in terms of six spectral densities. 

We finally note that in deriving the spectral representations (30) no use was made of the KMS condition. In this 
form, (with the spectral densities defined by (16), (26) and (28)), Eqs. (30) remain valid out of thermal equilibrium. 

IV. 4-POINT VERTEX FUNCTIONS IN QED IN THE HTL APPROXIMATION 

In this section, we calculate the seven 4-point vertex functions (10) for QED in the HTL approximation. In the 
HTL approximation QED has two non-zero 4-point vertex functions, shown in Fig. [1] below: 





ft 



Fig. [1] : 4-point vertex in QED (a - upper graph) and its cross term (b - lower graph) . 



The other two 4-point vertex functions (four-photon vertex and four-fcrmion vertex) are zero in the HTL approximation 
[3] . We calculate the diagrams in Fig. [1] as follows: the real-time Green functions are written in the Keldysh formalism 
using (5) and (10). The photon propagator in the Feynman gauge is given by —ig^D^) where D{K) is defined in 
(6). The fermion propagator is defined in (7). For the diagram in Fig. [la] we obtain 



where 



G^ a \ Pl ,P 2 ,P 3 ,P 4 ) = - 



G$ a) (P u P 2 ,P 3 ,P 4 ) 



d 4 K 



d 4 K 
(2^ 



H^(K)[aia 2 a 3 a 4 + r\f 2 f 3 f 4 + /i 626364 
+7-1/26364 + rif 2 / 3 a 4 + r 1 f 2 f 3 f 4 \ , 

if*" (IT) [01026364 + r x f 2 f 3 f 4 + hr 2 f 3 f 4 
+01/26364 + air 2 f 3 a 4 + air 2 r 3 / 4 ] , 



d 4 K 



G$ a) (P u P 2 ,P 3 ,P 4 ) = 
G^(P U P 2 ,P 3 ,P 4 ) = 
G?±\P 1 ,P 2 ,P 3 ,P 4 ) = 
GZ^(Pi,P2,P3,Pi) = 



(27T) 

d 4 K 



(27T)' 



d 4 K 



d A K 



d 4 K 
(2tt)4 



A H^ V {K) [ai6 2 6 3 6 4 + r 1 f 2 f 3 f 4 + ha 2 f 3 f 4 

+rifi¥ 3 r 4 + 0162/364 + 0162^/4] , 
■H^(K) [01626364 + r x f 2 f 3 f 4 + f\a 2 d 3 f 4 



+^1/263^4 + nf 2 f 3 f 4 + 016263/4] , 

H^ U (K) [01626364 + rif 2 f 3 r4 + /162/364 + 7-164/2/3 
+rif26 3 6 4 + f\d 2 f 3 f 4 + r 1 f 3 f 2 f 4 + a!a 2 f 3 f 4 } , 

if*" (if) [01626364 + rif 2 f3f 4 + 7-1626364 + 01^27=3^4 

+/1/26364 + hhr 2 a 4 + fif 4 f 2 f 3 ] , 
H^ v (K) [01626364 + r x f 2 f 3 f 4 + 0162637=4 + rif 2 f 3 a 4 
+/1/46263 + f 2 f 4 r x a 3 + hf 4 nf 2 ] , 



H^{K) = e 4 7 a(# - f 1 + m)^(f 2 + P 3 + $ + m) lv {P 2 + If + m) 7a 



(34a) 
(34b) 
(34c) 
(34d) 
(34c) 
(34f) 
(34g) 
(35) 



and we used the notation ai = Da(K), r\ = Dr(K), fi = Dp{K) for photon propagators and Oj = Da{Kj), 
fj = DniKj), fj = Dp(Kj) [j = 2,3,4 and K 2 =P 2 +K, K 3 =P 2 +P 3 +K, K 4 =K-P X \ for electron propagators. Note 
that the first two terms in each equation of (34) vanish by contour integration since all poles are located on the same 
side of the real axis. 

Before doing any momentum integrations, we can see immediately from (34) that the one- loop retarded functions 
{Cm, i — 1,2,3,4) are linear in the distribution functions, since they are linear in the symmetric propagators fi 
and/or fi. Terms with higher powers of the distribution functions are explicitly cancelled in the Keldysh formalism. 
This type of cancellation also occurs for the 2- and 3-point retarded- advanced functions [27]. In contrast, the mixed 
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retarded-advanced functions are bilinear in the distribution functions. One of the advantages of the Keldysh formalism 
is that it allows us to immediately identify the leading order temperature dependence. 

In the HTL approximation, the external momenta are soft (of order eT with e <C 1) and the loop momenta are 
hard (of order T) [4,5]. In this approximation, we can neglect the external momenta and the bare electron mass m 
relative to the loop momenta. For ^[^(K) one then obtains 

H^{K)^%e A K^K v If . (36) 

We first consider G"^. Fig. [la] gives 

GZ [a) {Pi, P 2 , P 3 , Pi) = \J 0^ H^iK^h^asai + nf 2 a 3 a 4 + nfa/ao* + TWa/i] 

= h+h+h+ h- (37) 



From (6), (7) we have 



(H(K) = ai{K) = — — -jl 



K 2 — mf — isgn(fc )e 
1 



n {K) = n(K) = — y— — - . 

K 2 - mf + isgn(k )e 

h{K) = -27nS(K 2 )(l + 2n B (k)), 

fi(K) = -2 m S(K 2 - m 2 )(l - 2n F (k)) , (38) 

where i—2, 3, 4 and m2=m 3 =mi=m which stands for the electron mass. We calculate the contribution from each of the 
terms in (37) in the HTL approximation: K ~ T 3> Pi, P2, P3, P4, m. For the temperature-dependent contributions 
we obtain 

h = J dkkn B (k) J dnV^V v Y ■ KsKK . (39a) 

I 2 = ^3 J dkkn F (k) J dnV»V"Y • 62636^, (39b) 
h = -^y J dkkn F (k) J dnV^VY • 63 64 633 > (39c) 
7 4 = --^3 J dkkn F {k) J dQ.V»V v Y -btbtbtz, (39d) 



where 

and V — (1, V) with V — k/k a light-like unit vector, and J dft denotes integration over the orientation of V. 
By interchanging P3 and P4 we obtain the contribution from Fig. [lb] to Gr\: 



with 



G Ri b) =Ji + J2 + J 3 + Ji, (41) 



Ji = 7^3- Jdkkn B (k) J dQV^V'y -b+b^b^, (42a) 
J 2 = ^3 y dkkn F {k) J dnV»V v Y -KK h t2 , ( 42b ) 



• 4 

le 



J 3 = 7^3- y dkkn F {k) J dnv»v v Y ■ ^KKa . (42c) 

J 4 = --^^ y dkkn F (k) J (KIVVY -b+b^b+ 2 . (42d) 
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Combining (39) and (42), the total contribution to from Fig. [1] reads 



g^ 1 (p 1 ,p 2 ,p 3 , p 4 ) = GZ {a) + gZ w = % ^f- J dn vvy ■ b 23 b 2 M - b 2 ) (43) 

where 

e 2 f e 2 T 2 



dkk{n B (k) +n F (k)) = — — (44) 



is the usual expression for the square of the dynamically generated thermal electron mass. 

Proceeding similarly we can calculate the other three retarded 4-point functions in finite temperature QED in the 
HTL approximation. The final results, including both diagrams in Fig. [1], are 



G% 2 (P U P 2 ,P 3 , P t ) = f dfi V»Vyb+ 3 b+ A (bI - b+) , (45a) 

' 2 2 

G^(P U P 2 ,P 3 , Pi) = J dQ VVYb&b&iK - b 2 ) , (45b) 

G% 4 (P U P 2 ,P 3 , Pi) = j dQ V»V v Yb 23 b 2i (b^ - b 2 ) . (45c) 

Our results agree with those obtained in the ITF [18]; to leading order, all of the retarded 4-point functions are 
proportional to T 2 . 

Next we consider the mixed retarded-advanced 4-point functions. We begin with G^ ra . The leading order 
temperature-dependent piece from Fig. [la] is 

G^a\Pu P2, P 3 , Pi) = \f 0ji H^{K)[ha 2 h<ii + nhhai + /ia 2 f 3 / 4 + nf 3/2/4] 

= T 1 +T 2 +T 3 + T i , (46) 

where the functions Ti ; i = {1,2,3,4} represent the contributions of the four terms respectively. Inserting the 
propagators (38) and using the HTL approximation we find 



T 2 
T 3 
T 4 

where 



(2tt) 3 

e 4 



cp J dQ V^Yb+b 2 (b 23 - b+ 3 ) , (47a) 

ad„ J dnv^Yb+^ib, b+) , (47b) 

j^c fj j dCl V»V»Yb 2 bt 3 (bt for) , (47c) 

j^d J dnv^Ybtb 2 (bt 2 b u ) , (47d) 



cp = J dkk[n B (k)(l—n F (k)) - n F (k)(l+n B (k))} , 

d = -2 J dkkn F {k){l-n F (k)). (48) 

Inserting the standard forms for the Bose-Einstein and Fermi-Dirac distributions we obtain in thermal equilibrium 

C(8 =0, dp = -(21n2)T 2 . (49) 
The contribution from Fig. [lb] to G^ ra is 

Grara = S 1 + S 2 + S 3 + Si (50) 

with 
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Si = 


g 4 

(2tt)3 °n 


1 dnv^Vfb+b^ib^ 




(51a) 


s 2 = 


e 4 

J 




-bt), 


(51b) 


s 3 = 


e 4 


1 dVLV^V u fb^b+ A {b+ 




(51c) 


s 4 = 


(27r) 3 d/3 J 


f dn v^Vfb+bzib+z 


K 2 ) ■ 


(51d) 



Adding (47) and (51) we obtain 



^ vara 



(27T) 



3 / dQ V»V»Y dp 6 2 -(&+(& 3 - + bl) - b^(b+ + &+)) + c p 6 2 -(6+(6 2 - 3 + fe 2 - 4 ) - bi(b+ + &+)) . (52) 



In spite of (49) we keep the terms ~ cp because we will discuss the extension of these results to non-equilibrium 
situations in Sec. VII. 

The other two mixed retarded-advanced 4-point functions are similarly calculated: 



QfXV 

raar 



rraa 



(2tt) 3 

e 4 
(2tt) 3 



dn WV V Y dp b+(bt 2 (bz + 6 4 ) - b^(b+ + b+)) + cp 6 2 -(&+(6 2 - 3 + 6 2 - 4 ) - + b 23 )) 



dQ V»V V Y cp (6 2 -6+(6 2 -3 + 6 2 - 4 ) - b+bi(b+ 4 + 6+)) 



(53) 



The last function vanishes in equilibrium since cp = 0. The other mixed retarded-advanced 4-point functions are 
again proportional to T 2 . Furthermore, since V 2 =Q, all seven 4-point functions are traceless (Gfr = 0) in the HTL 
approximation. 

V. 4-POINT SPECTRAL FUNCTIONS IN HOT QED IN THE HTL APPROXIMATION 

From (43), (45), (52) and (53) one obtains the following spectral representations for the 4-point retarded-advanced 
functions in the HTL approximation: 



G^(W1,W2,W3,W 4 ) 



(2^ 



r 2 <*>2, UI3 



dnidfl2dn 3 dfl4 (5(^1 + ^2 + ^3+^4) Phtl^ 1 ' ^3' a 23 a 2i{ a \ ~ a 2 ) ' 
,W 4 ) = y'^1^2^3^4^(^l+^2+^3+^4)p^ L (^l,^2,^3,^4)a23 a 24( a r _ a t) ) 

w 2 , w 3 ,w 4 ) = —^-3 J dn 1 dn 2 dn 3 dn 4 <5(f2i+n2+0 3 +Q 4 ) /9^ L (Oi,n 2 , 3 , 4 ) a^a^aj" - a 2~); 
G^(wi,W2,w 3 ,a>4) = ^ 2 ^ 3 y ^1^2^3^4(5(^1+^2+^3+^4) (Oi , 2 , 3 , 4 ) a^a^af -a^), 

dOidf2 2 ^3^4 (5(01 + 02 + ^3 + ^4) Phtl(°1' Q 2, 3 , 4 ) 
*1 a 2 ( a 12( a 3 + a 4 ) - a r2( a 3 + a 4 )) + a 2 («2 (d^ + O^) - a 2 ai~( a 24 + a 23)) 
G^ ar (wi,W2,W 3 ,W4) = ^3 y dr»id02d0 3 d0 4 J(Oi + 02 + 03 + 04)^ L (Oi,02,0 3 ,r! 4 ) 

• "lOi" (0^(03 + a 4 ) - af 2 (<4 + aj)) + o^af (a^ + a 24) - a 2~ a r( a 24 + a 23) 

G^ oa (wi,w2,w3,w 4 ) = 7^ / rfr>ido 2 do 3 do 4 5(Oi+o 2 +o 3 +r! 4 )p^ L (r!i,o 2 ,o 3 ,o 4 ) 



G rara( W l> W 2,W3,W4) = 



(27T) 



(2tt) 3 ., 

a 2 (a 2 (a^ + a 24) - a2~ a i~( a 24 + a 23)) 



(54a) 
(54b) 
(54c) 
(54d) 

(54c) 
(54f) 
(54g) 



where 
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^(n 1 ,n 2 ,n 3 ,n 4 ) = 2w 2 e 2 m} J dfiv^YS(P 1 -v)S(P 2 -v)S(p 3 -v) (55) 



(with Pf = Sli, P t = P t ; i={l, 2, 3} on the r.h.s.) and 

= dpe 4 _ = c^e 4 
2iTT 2 e 2 m 2 i ' 2iTT 2 e 2 m 2 j ' 

These equations also hold in the non-equilibium situations studied in Sec. VII. In equilibrium we use (49) to obtain 

81n2 , N 

ai = i — 5— ; a 2 = (57) 

Note that the frequencies Oj should not be confused with the angular factor 0. From V 2 =0 it follows that 

(/W)£ = 0. (58) 

In contrast to the general spectral representations (30) of the 4-point functions given in Sec. Ill, in the HTL 
approximation the spectral representations (54) involve only a single, real- valued spectral density p^ L - This agrees 
with the result of Taylor [18] who showed for QCD in the ITF that in the HTL approximation all spectral densities 
degenerate to a single function. Clearly the spectral density (55) is symmetric in the Lorentz indices and under 
interchange of the external momenta. Using the 5-functions in (55) it is also easy to show that is transverse with 
respect to all external momenta: 

^p£l=°. i = 1,2,3,4. (59) 
A previous RTF calculation of the 3-point spectral function in hot QED in the HTL approximation gave [26] 

p£ TL (fii,fi 2 ,fi 3 ) = — 2^ / dSlWY S(P 1 -V)6(P 2 -V). (60) 
By comparing (55) with (59) we find the "sum rule" 

dn 3 pf(Q u n 2 , o 3 , ~(n 1 +n 2 +n 3 )) = 4^ e p/tisiuth, -(Oi+o 2 )) . (6i) 



This agrees with the ITF results in [18]. A similar sum rule exists between the 2- and 3-point spectral densities: the 
2-point spectral density is given by 

p HTL (fii,-fii) = ^ J dQYm-V), (62) 
and the temporal component of the 3-point spectral density p^ TL O2, —(^1+^2)) obeys the sum rule 

dn 2 p HT ^n 1 ,n 2 ,-(Q 1 +n 2 )) = 7rep HT jn 1 ,-n 1 ). (63) 

We believe that analogous sum rules exist for the higher order n-point functions. 



VI. WARD IDENTITIES BETWEEN 3- AND 4-POINT FUNCTIONS IN THE HTL APPROXIMATION 

In this section, we verify the HTL Ward identities between the 3- and 4-point vertex functions in QED, without 
doing momentum integrations. We split propagators in the integral expressions (34) by using the identities [28] 

D R/A {K 2 ) D R/A (K 3 ) = 1 (D R/A {K 2 ) D r/a {K 3 )) (64) 

This trick allows us to rewrite (34) as 
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G m =e I 7^4 op Ui\ a 2 - a3)«4 + rif 2 a i - ri/ 3 a 4 + n(r 2 - r 3 )/ 4 J 



^(°) _ „4 

7?9 — e 



(2tt) 4 2P 3 -if 
f d 4 K AK^K U % 



f d"K AK^£ - - 

J (27T) 7 2P 3 -if ~ r3 ^ 4 ai -^ 2Ct4 ~ ai /3«4 + ai(r 2 - r 3 )/ 4 ] , 

^(») 4 [ d 4 K 4K»K»If - -_ - 



G' 



„ (o) _ 4 [ d^AK^i£Ll_ 

) 4 2P 3 • A 



«4 



y (2tt 



nnv(a) _ 4 



d 4 if AK^K v Jf 



[/i(a 2 - a 3 )r 4 + ri/ 2 r 4 - ri/ 3 r 4 + oi(a 2 - a 3 )/ 4 ] , 



-hha-i + fi(r 2 - a 3 )a 4 + /i(a 2 - r 3 )/ 4 + ri/ 2 / 4 + ai(a 2 - r 3 )r 4 ] , 



(2tt) 4 2P 3 • K 

G%£> = e 4 | ^ A 2P^K [ri{R2 ~ S3) " 4 + ai(f2 ~ f3)r " 4 + /l/2S4 " /l/3 " 4 + /l/4(f2 " f3)] 

/(Pif AK^ K v If 
~{2nj 4 ~ 2P 3 -K ^ l( -° 2 ~ Ss ' )f4 + ri( - f2 _ f3 ^ 4 + ^ 4 ^ 2 ~ + fof iri ~ hhr\\ 

where we have used the HTL approximation to write K\ — K\ = 2K ■ P 3 + P| w 2_ftT • P 3 . 



(65a) 
(65b) 
(65c) 
(65d) 
(65e) 
(65f) 
(65g) 




K-Pl 



P2 



Fig. [2]: 3-point vertex in QED 

The integral expressions for the 3-point vertex functions in QED are obtained from Fig. [2]. In the Feynman gauge 
we have 



G^(P 1 ,P 2 ,P 3 +P 4 ) -e S 



G^(P l5 P 2 ,P 3 



+P4) = e 3 | 



d A K 
(2^ 
d 4 K 

(27T) 4 

d 4 K 



(2i4r AI ^)[aia 2 a4+rir 2 r4+/ir 2 r4+ai/ 2 a4+air 2 / 4 ], 
(2Ar A1 ^)[aia 2 a 4 +?'ir : 2 f4+/ia 2 a4+ri/ 2 a4+rir 2 /4], 



G^ (Pi,P 2 ,P 3 +P 4 ) -e 3 / 7 ^T(2^^)[aia 2 a4+r 1 f 2 F4+/i0 2 f4+r 1 / 2 f4+a 1 a 2 / 4 ], 



G^(Pi,P 2 ,P 3 +P 4 ) = e E 
G£ o 0Pi,P 2 ,P 3 +P 4 ) = e E 



(2^) 4 
d 4 £T 



(27T) 
rf 4 if 

(2^ 
d 4 K 



4 (2K fl If )[a 1 a 2 r 4 +r 1 r 2 a4+r 1 f 2 f 4+ f ^2/4], 
(2if At ^)[aif 2 a 4 -l-ria 2 f4-l-ai/ 2 /4+/i/ 2 r4], 



(27T) 



4 (2K tl If) [ria 2 a4+air 2 r4+/i/ 2 a 4 +/ir 2 / 4 ] , 



(66a) 
(66b) 
(66c) 
(66d) 
(66e) 
(66f) 



where for later convenience we have chosen to write on the l.h.s. — (Pi+P 2 ) = P3+P4. Comparing (65) and (66) we 
obtain the following Ward identities between the vertices shown in Fig. [la] and Fig. [2]: 
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P 3 „G R \ {a) {P u P 2 , P 3 , P t ) = e(G Ri (P u P 2 , P 3 +P 4 ) - G v m {P l ,P 2 +P 5 ,P i )) , (67a) 

P^G R v 2 {a) {P u P2, P 3 , Pi) = e(G R (P 1 ,P 2 ,P 3 +P 4 ) - G» R (P 1 ,P 2 +P 3 , P 4 )) , (67b) 

P 3 „G R v 3 {a) (P u P 2 ,P 3 , P 4 ) = e(G^ (P 1 ,P 2 ,P 3 +P4) - G R (P 4 ,P 2 +P 3 ,P 4 )) , (67c) 

P 3tl G R v £ a) {P u P 2 , P 3 , P 4 ) - e(G^ (P 1 ,P 2 ,P3+P 4 ) - G Ro (P 1 ,P 2 +P 3: P 4 )) , (67d) 

P 3 „G?±\P U P 2 ,P 3 , P 4 ) = e (G^(P 1 ,P 2 ,P 3 +P 4 ) - G^ (P, P 2 +P 3 , P 4 )) , (67e) 

^V^± } (Pi, P 2 , P 3 , P 4 ) = e(G v F {Pi,P2,P3+Pi) - G F (P U P 2 +P 3 ,P 4 )) , (67f) 

P Sfl G?±\P u P 2 ,P 3 , P 4 ) - e(G^ (P 1 ,P 2 , P 3 +P 4 ) - GJ- (Pi,P 2 +i%,P4)) • (67g) 
Similarly one can obtain the Ward identities satisfied by the diagrams in Fig. [lb] and Fig. [2]: 

P^G^iP^P^P^P,) = e(G Rl (P 1 +P 3 ,P 2 ,P 4 ) - G Rl (P u P 2 ,P 3 +P 4 )) , (68a) 

P^G R ^ h) {P u P 2 , P 3 , P 4 ) = e(G£(Pi+i^,P 4 ) - G R (P U P 2 , P 3 +P 4 )) , (68b) 

P^fHP^P^P^P,) = e(G Ro {P x +P 3 ,P 2 , P 4 ) - GUA,^2,P3+P4)) , (68c) 

P^G R f ] (P l: P 2: P 3 ,P 4 ) = e(G Rl (P 1 +P 3 ,P 2 ,P 4 ) - G Ro (P u P 2 ,P 3 +P 4 )) , (68d) 

P3 M G^W(P 1; P 2 ,P 3 ,P 4 ) = e (G^(Px+P 3 , P 2 , P 4 ) - G F {Pi,P2,P3+Pi)) , (68e) 

^G^W (P 4 ,P 2 ,P 3 ,P 4 ) = -eG F (P 1: P 2 ,P 3 +P 4 ) , (68f) 

P3 M G^W(P 1; P 2 ,P 3 ,P 4 ) = e^JPi+i^Pa.ft) - G£ (P 1; P 2 ,P3+P 4 )) . (68g) 

By combining (67) and (68) one obtains the following HTL Ward identities between the 3- and 4-point vertex functions 
in QED in the RTF: 

P 3 ^G R \(P U P 2 ,P 3 ,P 4 ) = e (G^(P+P 3 , P 2 , P 4 ) - G Ri (P u P 2 +P 3 ,P 4 )) , (69a) 

P 3 ^ 2 (P U P 2 ,P 3 ,P 4 ) = e(G R (P 1 +P 3 ,P 2 ,P 4 ) - G R (P U P 2 +P 3 ,P 4 )) , (69b) 

P 3 ^ 3 (P U P 2 ,P 3 ,P 4 ) = e(G R (P 1 +P 3 ,P 2 ,P 4 ) - G Ro (P 1 ,P 2 +P 3 ,P 4 )) , (69c) 

P^G R \(P 1: P 2 ,P 3 ,P 4 ) = e(G Ro (P 1 +P 3 ,P 2 ,P 4 ) - G Ri (P u P 2 +P 3 ,P 4 )) , (69d) 

P 3 »GZ-a(PuP2,P3,P4) = e(G F (P 1 +P 3 ,P 2 ,P 4 ) - G Fo (P u P 2 +P 3 ,P 4 )) , (69e) 

P 3 ^ ar (P u P 2 ,P 3 ,P 4 ) = -eG F (P u P 2 +P 3 , P 4 ) , (69f) 

P 3 ^ aa (P u P 2 ,P 3 ,P 4 ) = e(GUA+^3, P2, P 4 ) - G Fo (P l ,P 2 +P 3 ,P 4 )) . (69g) 



Notice that in equilibrium the last equation gives simply = (see (49) and (53)). These results are structurally 
similar to the zero temperature Ward identity [29] and to the Ward identities between the retarded 3- and 4-point 
HTL vertices in QCD obtained within the imaginary time formalism [8] or from kinetic theory [10]. 

VII. EXTENSION TO NON-EQUILIBRIUM SITUATIONS 

The HTL method has been widely used in qualitative studies of the QGP and in other explicit calculations at finite 
temperature [5]. This powerful method was derived within the ITF for equilibrium field theory, and this is where its 
range of applicability can be clearly defined [4,5,8]: HTL resummation is required and applicable for the study of soft 
processes in a weakly interacting plasma, with momentum scale p~gT<T(g<l) much below the typical thermal 
or "hard" momenta in the medium. Realistic physical systems, on the other hand, arc frequently out of equilibrium, 
which means that calculations must be carried out in the RTF. As long as the local momentum distribution f(x,p) 



13 



doesn't deviate too strongly from a thermal one, and can still be roughly characterized by a parameter T which gives 
the scale for the typical momenta at point x, a similar "hard loop" (HL) resummation scheme should be applicable for 
the study of the dynamics of soft modes (p~jT<T) in a weakly interacting (j « 1) non-equilibrium system [28]. 
This approach is expected to be reliable for equilibration processes requiring momentum exchanges which happen 
slowly, on a time scale t ~ (t^T) -1 » l/p°. This reasoning does not apply to arbitrary non-equilibrium situations, 
but only to plasma states which are sufficiently close to global thermal equilibrium; on the other hand, deviations 
from chemical equilibrium can be arbitrarily strong without invalidating the scheme. 

In this section we discuss the generalization of our results to this particular type of non-equilibrium scenarios. 
Since we have used the Keldysh or closed time path contour, the procedure is straightforward. Equations (1) and 
(2) can be used in non-equilibrium situations by replacing the x-independent equilibrium Bose-Einstein and Fermi- 
Dirac distributions by non-equilibrium distributions (Wigner functions) and Jf(x,p) which depend on the 

space-time coordinate and the four momenta [21]. Equations (6) and (7) remain valid for bare propagators, with the 
equilibrium distribution functions (ns^np) replaced by the non-equilibrium ones (/b,/f), but the last equations in 
each of (6) and (7), (i.e. the fluctuation-dissipation theorem which reflects the KMS condition), no longer hold for 
full propagators [28]. 

Equations (30) give the spectral representations of the seven 4-point functions (10) in equilibrium. These expressions 
remain valid out of equilibrium since the KMS condition was not used in their derivation. Out of equilibrium there 
are eight more 4-point vertex functions {G arrri G rarr , G rrar , G rrrai G ara r, G aarr , G arra and G rrrr ) with similar 
spectral representations. They can be derived using the same methods as those given in Sec. III. 

The results in Sees. IV and V can be generalized to non-equilibrium situations by simply replacing in (44) and (48) 
the Bosc-Einstcin and Fcrmi-Dirac distributions nsip ) and npip ) by non-equilibrium Wigner functions Jb{x,p) and 
/f (x,p). The thermal mass mp in (44) and the factors cp, dp in (48) become functions of x, m(x), c(x), d(x), and it is 
only in the special case of thermal equilibrium that c = as obtained in (49). These generalized results correspond to 
the "hard loop" (HL) approximation in slightly off-equilibrium plasmas and are applicable in the situations discussed 
at the beginning of this section. We have shown that the 7 components of the 4-point function given in (10) correspond 
to 5 spectral densities which degenerate in the HL approximation to a single real function. This result remains valid 
out of equilibrium, where the KMS conditions do not hold. The eight additional 4-point functions which must be 
included out of equilibrium will involve more spectral functions, but in the HL approximation it is straightforward 
to show that they degenerate to the same HTL spectral density (55), with the thermal mass mp replaced by the 
generalized mass m(x). 

The one-loop Ward identities in Sec. VI were obtained without explicitly doing any momentum integrals; the result 
thus holds independently of the form of the distribution functions (which enter through the functions fi, fi in (65) and 
(66)). Out of equilibrium there will be additional Ward identities for the eight additional vertex functions which can be 
derived in the same way. For example, the Ward identities that correspond to the three functions G aarrj G arra7 G arar 
have been calculated, and they have the same form as those obtained in (69 e,f,g). 

VIII. CONCLUSIONS 

We have studied the 4-point vertex functions in the Keldysh formulation of real-time thermal field theory [22] . This 
formalism has recently gained increased popularity because it avoids the need for analytical continuation associated 
with the imaginary time formalism, and it allows for a generalization to non-equilibrium situations. Starting from the 
largest time and smallest time equations which relate the 16 components of the 4-point real-time vertex function, we 
have derived spectral integral representations for the 7 retarded-advanced 4-point functions (10). We have explicitly 
calculated these 7 functions for hot QED at 1-loop order in the HTL approximation, using the Keldysh representation. 
In particular, we have demonstrated the usefulness of this representation for the calculation of many-point functions 
and showed how terms with high powers of the distribution functions cancel explicitly before any momentum inte- 
grations are done. The fully retarded vertex functions are found to be linear in the distribution functions [27]. In the 
HTL approximation, they have the same structure as their QCD counterparts which were derived in the ITF [4,30] as 
well as from kinetic theory [10]. We have also calculated the three mixed retarded-advanced 4-point functions which 
are needed for a complete description of the real-time 4-point vertex and which have not previously appeared in the 
literature. They are quadratic in the distribution functions. One of them (G rraa ) was found to vanish in thermal 
equilibrium, while the other 4-point functions are all proportional to T 2 . 

We have derived spectral integral representations for the 7 retarded-advanced 4-point functions (10). The spectral 
densities were calculated explicitly for hot QED in the HTL approximation and found in this case to degenerate to 
a single real spectral function. Its trace over the Lorentz indices vanishes, and it is transverse with respect to all 4 
external momenta. It obeys a simple sum rule with the correspondingly approximated 3-point spectral density. 
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By contracting the 7 retarded-advanced 4-point functions (10) for QED in the HTL approximation with one of the 
photon momenta we have derived one-loop Ward identities. This calculation was done without doing any momentum 
integrations, by comparing the integrands of the contracted HTL 4-point vertices with those of the HTL 3-point 
vertices. The resulting set of real-time Ward identities at finite temperature generalizes the zero temperature Ward 
identity and can be compared with the HTL Ward identity from the ITF and from kinetic theory. Due to the matrix 
form of the real-time thermal Green functions, the Ward identities have a more complex structure: one finds a whole 
class of finite temperature Ward identities which relate retarded-advanced 4-point functions to retarded-advanced 3- 
point vertex functions. If either the ingoing or the outgoing fermion leg has the largest time, the Ward identity involves 
only one kind of retarded-advanced 3-point vertex functions, and the right hand side vanishes when the contracted 
photon momentum (P 3 ) goes to zero. If one of the photon legs has the largest time, two different retarded-advanced 
3-point functions are involved, and the right hand side no longer vanishes (Eqs. (69c-f); the corresponding component 
of the 4-point vertex is thus singular in the zero-momentum limit P3 — > [31]. Similar features were found for the 
Ward identities relating the HTL 3- and 2-point functions in hot QED [20] . 

Our Ward identities are not general, in the sense that they were not directly derived by functional methods from 
the RTF path integral. However, in the case of the 3-point functions it is known that the one-HTL Ward identities 
[20] are identical in structure with those satisfied by the exact vertex functions [32] . We therefore expect the same to 
be true in the present case. 

Finally, we have discussed the generalization of our results to certain classes of "moderately non-equilibrium" 
situations (Sec. VII) in which a generalized HTL approximation ("hard loop approximation") can be applied. Out of 
equilibrium the number of independent 4-point functions is larger, but in the HL approximation the spectral functions 
have the same structure as in equilibrium. All of the nice features of equilibrium HTLs persist in such "moderately 
non-equilibrium" situations. 
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